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Abstract 

We review the accelerating (mainly, dark energy) cosmologies in modified grav- 
ity. Special attention is paid to cosmologies leading to finite-time future singulari- 
ties in F{R), F{G) and J^{R, G) modified gravities. The removal of the finite-time 
future singularities via addition of i?^-term which simultaneously unifies the early- 
time inflation with late-time acceleration is also briefly mentioned. Accelerating 
cosmology including the scenario unifying inflation with dark energy is considered 
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in F{R) gravity with Lagrange multipliers. In addition, we examine domain wall 
solutions in F{R) gravity. Furthermore, covariant higher derivative gravity with 
scalar projectors is explored. 

PACS numbers: 04.50.Kd, 95.36.+X, 98.80.-k 
1. Introduction 

It is observationally implied that the current expansion of the universe is 
accelerating. Provided that the universe is homogeneous, two representa- 
tive approaches to account for the current cosmic acceleration exist. The 
first is to assume the existence of the so-called dark energy whose pressure 
is negative (for a recent review, see, e.g., P^). The second is to consider 
that a gravitational theory would be modified at the large distance scale. 
The simplest theory is F{R) gravity (for reviews, see, for example, [2])- 

In this paper, we examine the accelerating (dark energy) solutions of 
modified gravity which may produce future singularities. We concentrate on 
reviewing the results in Refs. [3l HI 13 [6] on theoretical aspects of modified 
gravity theories with presenting dark energy components. In particular, 
we study the finite-time future singularities in F{R), F{G) and T{R,G) 
gravity theories O d [8], where R is the Ricci scalar, G = R^ - AR^yRf^"" + 
Rfj.upaR'^'^'"^ with R^i, and R^upa the Ricci tensor and the Riemann tensor, 
respectively, is the Gauss-Bonnet invariant, and J-^{R, G) is an arbitrary 
function of R and G. This is a generalized gravity theory including both 
F{R) and F{G) gravity theories. We also discuss the removal of the finite- 
time future singularities in F{R) gravity via addition of i?^-term which 
simultaneously leads to the unification of early-time inflation with late- 
time acceleration [9j. In the frameworks of F(G) or J^(R,G) theory, the 
corresponding term may be different, of cours e We note that as related 
studies, the finite-time future singularities [IO||ll((12J and the realization of 
the phantom phase including the crossing of the phantom divide [13j have 
also been examined. Furthermore, the features of the finite-time future 
singularities in non-local gravity [H], modified teleparallel gravity [15] and 
its extended analysis in loop quantum cosmology (LQC) [16] has recently 
been investigated. In addition, dark energy in the context of F{R) gravity 
with Lagrange multipliers^^] is considered. We also present domain wall 
solutions in F{R) gravity [^7 Moreover, covariant higher derivative gravity 
with scalar projectors [6J is explained. We use units oi = c = h = 1 
and denote the gravitational constant SvtGn by = Svr/Mpi^ = 1 with 
the Planck mass of Mpi = G^^^^ = 1.2 x lO^^GeV. 

The paper is organized as follows. In Section 2, we explore accelerating 
cosmologies leading to the finite-time future singularities in F{R), F{G) 
and F{R, G) gravity theories. In Section 3, we study dark energy in the 
framework of F[R) gravity with Lagrange multipliers. In Section 4, we 
examine domain wall solutions in F{R) gravity. In Section 5, we inves- 
tigate covariant higher derivative gravity with scalar projectors. Finally, 
conclusions are presented in Section 6. 
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2. Finite-time future singularities in F{R), F{G) and J-'{R, G) 
gravity theories 

2.1. T{R,G) gravity 

The action of I'{R,G) gravity is S = J df^x^f^ (2^2) + Cm], 

where g is the determinant of the metric tensor g^y and Cm is the matter 
Lagrangian. This is a generic theory including both F{R) and F{G) gravi- 
ties. We take the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) met- 
ric ds"^ = —dt^ + a?{t) X]i=i 2 3 (f^^;*) . The Hubble parameter is given by 
H = a /a, where the dot denotes the time derivative of d/dt. In the FLRW 
background, with the gravitational field equations we find that the effective 
(i.e., total) energy density and pressure of the universe read peflf = Sk'^H'^ 

and Peff = —K'"^ (2H + 3H^], respectively. For the action in Eq. ([T|, we 



obtain 

\iFRR -F)- 6HTr + GF,G - 2iH^TG } , (1) 
(F^rR -F) + AHTr + 2Tr - GF,G 









PcS = 






-feflf = 











+ 16H{H + H'] Tg + SH'Tg 



} , (2) 



where F^r = dF{R, G)dR and F^g = dF{R, G)dG, and pM and Pm are 
the energy density and pressure of matter (which has been assumed to be 
a perfect fluid). 

2.2. Finite-time future singularities 

Provided that the Hubble parameter is written as 

^=(^7^ + ^- (3) 

where hs{> 0), ts{> 0), Hs{> 0), and 0) are constants, tg is the time 
when a finite-time future singularity occurs, and < t < ts- In what 
follows, we consider the case of Hg = 0. We note that even if /3 < and 
/3 is a non-integer value, in the limit t ^ ts some derivative of H diverges 
and hence the scalar curvature becomes infinity [lOj. Moreover, since the 
case of /3 = leads to a de Sitter space, we suppose /? 7^ 0. 

The finite-time future singularities are classified into four types [17] . 
Type I ("Big Rip" [IS]): In the limit i — > tg, a — ^ 00, peff 00 and 
|Peff| — ^ 00. The case that peS and Pes are finite values at t = is [H] is 
included. This happens for (3 = 1 and /3 > 1. In this paper, we regard 
the singularities for /3 = 1 as "Big Rip" and those for /3 > 1 as "Type I", 
(ii) Type II ("sudden" [20]): In the limit t ^ tg, a ^ Ug, PeS — ^ Ps and 
|Peff| ^ 00. This occurs for — 1 < /3 < 0. (iii) Type III: In the limit t ^ tg, 
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a — )• Os, Peff — ^ oo and |-Peff| oo. This appears for < /3 < 1. (iv) Type 
IV: In the limit t ^ tg, a ^ Og, PeS — ^ 0, |-Peff| — ^ 0, and higher derivatives 
of H diverge. The case that p^s and/or |Pcfr| become finite values at t = tg 
is also included. This is realized if /3 < —1 but /3 is not any integer number. 
Here, 03(7^ 0) and ps are constants. 

2.3. F{R) gravity with finite-time future singularities 

By taking F{R,G) = F{R), the action in Eq. ^ becomes that of F{R) 
gravity. With the method to reconstruct modified gravity I^H [22] , for 
the Hubble parameter to be represented in Eq. ([3|) we explore F{R) gravity 
models in which finite-time future singularities can appear. By introducing 
two proper functions -P(0) and (5(0) of a scalar field 0, which we regard 
as the cosmic time t, we rewrite the term J-{R, G) = P{t)R + Q{t) in the 
action in Eq. ([T|. In this case, by varying the action with respect to t we 
acquire {dF{t)/dt) R + dQ{t)/dt = 0. In principle, by solving this equation 
we have the relation t = t{R). If we substitute it into the above form of 
F{R,G) = P{t)R + Q{t), we find F{R) = P{t = t{R))R + Q{t = t{R)) 
and hence the original action is found again. We express the scale factor as 
a{t) = aexp {g{t)) with a a constant and g{t) a proper function. Here, we 
neglect the contribution from matter because when the finite-time future 
singularities appears, the energy density of dark energy components are 
completely dominant over that of matter. In this case, the gravitational 
field equations yield 

p{t) - mm + m)p{t) = o , Q{t) = -e [m)? pit) + mpit) ■ 

(4) 

Accordingly, if we find the solutions P{t) and Q{t) of these equations, by 
plugging those into F{R) = P{t)R + Q{t) with t = t{R) we obtain the 
concrete form of F{R). We acquire the followings consequences. 

(a) /3 = 1 [Big Rip]: For /is > 5 + 2^/6 or /i^ < 5 - 2^/6, F{R) oc i?« with 
q = (1/4) (3 + /is + V^s - 10/is + 1) , whereas if 5 - 2^6 < K < 5 + 2a/6, 



F{R) oc i^C^^+i)/^ X (Oscillating part). 

(b) /3 > 1 [Type I]: F{R) oc exp |(/is/ [2 (/3 - 1))] (^)^''"'^^^'''^| i^-V^x 



(Oscillating part). 



.(/3-l)/(/3+l) 



(c) < /3 < 1 [Type III]: F{R) ~ exp (-6/3/isi?)^ 

(d) /? < [Type II (-1 < /3 < 0) and Type IV (/3 < -1 but /3 is not an inte- 
ger)]: F{R) ~ (-6/is/3i?)(^'+2^+^)/[«(^+^)l exp {-6hm^^''^^^^^'^' 
Here, "~" means the asymptotic behavior in the limit t — )• is- 

It is remarkable that adding i?^-term to such a theory, one removes fu- 
ture singularity. (Note that R^ gravity was proposed as inflationary model 
in Ref. p3] (celebrated Starobinsky inflation) and was used for the first 
unified inflation-dark energy modified gravity proposed in Ref. [9j). Hence, 
we not only remove singularities by adding R^ term but also unify the dark 
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energy era with inflation in such a way (for realistic models of such unifi- 
cation, see [2]). Several viable F{R) gravity models which unify inflation 
with dark energy and do not contain the finite-time future singularities are 
listed below [2j: 

FiR) = R+ ''fl-'f +b,R\ (5) 



F{R) = - 2A 1 1 - e-^/^"^^^ I - Ai 1 1 - e-W«iJ' I + b^RT^^-^^R- , (6) 

with hj {j = 1,...,4), 65(> 0), &6(> 0) and / constants. In Eq. r(> 

1) is a natural number, and Ri are constants, Ri and Ai are transition 
curvature and expected cosmological constant at the inflationary stage, 
respectively 



2.4. F{G) gravity w^ith finite-time future singularities 

With the same method as in F(R) gravity in Section 2.3, it is possible to 
execute the reconstruction of F{G) gravity models in which the finite-time 
future singularities occur. The action of F{G) gravity [25] is described by 
Eq. ([T]) with F{R,G) = R + F{G). In this case, the gravitational field 
equations yield 

2^ (p\t)P{t)) - 2l\t)P{t) + m = , Q{t) = -2Al\t)P{t) - 6-g\t) . 
The results are as follows [3]. 

(a) /3 = 1 [Big Rip]: For / 1, F{G) = [^6hl{l + h,)/ %{! - /is)]} VG+ 

c^G{/is+i)/4+c2Gwithci andci constants. If /is = 1, F{G) = ^VGln{-fG) 
with 7(> 0) a positive constant. 

(b) /3 > 1 [Type I]: F{G) = -VgVG. 

(c) < /3 < 1 [Type III], -1/3 < ^ < [Type II], -1 < /3 < -1/3 [Type 
II] and P < -1 (but /3 is not integer) [Type IV]: F{G) = 6/1^(3/3 + + 
l)-i [|G|/(24/.3|/?|)]2^/(3^+i). 

(d) /3 = —1/3 [Type II] (this is a special value in this case): F{G) ~ 

[1/ (4\/6/i3)] G{G + Sh^y^^ + (2/V6) {G + Sh^^^^^ We remark that the 
finite-time future singularities appearing in the limit G — )■ ±oo can be 
removed by the additional term diG^, where di{^ 0) is a constant, and g > 
1/2 and Q ^ 1. Furthermore, the finite-time future singularities emerging 
in the limit G — ?■ 0~ can be cured by adding the term diG^, where q{< 0) 
is an integer [3]. 

2.5. F{R,G) gravity with finite-time future singularities 

Using the similar procedure in Section 2.3, we reconstruct the form of 
F{R, G) leading to the finite-time future singularities. With proper func- 
tions P{(p), and Q{(p) of a scalar field (p, which we indentify with t, we 
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represent the term F[R, G) in the action in Eq. ([T]) as P{t)R+ Z{t)G+Q{t). 
Varying this action with respect to i, we find {dP{t) / dt) R + {dZ{t)/dt) G + 
dQ{t)/dt = 0. By solving this equation, we obtain t = t{R,G). Combin- 
ing this and the above representation P{t)R + Z{t)G + Q{t), we acquire 
T{R,G) = P{t)R + Z{t)G + Q{t). It fohows from the gravitational field 
equations, the conservation law, a{t) = aexp (g(t)), and H(t) = g{t) that 

d'^Pi't) .-^^^d^Zit) .,,dP{t) 

+ 4 (2^1 - -9\t)) ^ + 2^(t)P(t) = , (8) 
Q{t) = -6 (4f (i)^ - i^Ht)P{t) - . (9) 

For P{t) / 0, T{R, G) can be described as T{R, G) = Rg{R, G) + /(i?, G) 

with g{R, G){^ 0) and f{R, G) generic functions of R and G. We show the 
results 

(a) ^ = 1 [Big Rip]: For < /is < 5-2^/6 or /is > 2+V6, we find T{R, G) = 
aiR'i+ + asfi"- + with q± = (1/4) (^3 + /is ± V^i - lO^s + l) • 
Here, ai, a2 and 5 are constants. There also exists the following model: 

where /(i? G) — / ~a(a:^+2)fi±Va^(a'+2)2R2+24fes[4fcgfa+a(3;+2+fes)](3:+4)<5zG \ 

Here, a and x are constants. 

(b) /? > 1 [Type I]: F{R,G) = -AhlXf{R,G)R + X{f{R,G)f+'^^G + 

1/(2/3) 

, where 



24/i4A(/(i?,G))i-2'^ with/(i?,G) 



(l/2+/3)G 



A is a constant. 

(c) /3 < 1 [Type II (-1 < /3 < 0), Type HI (0 < /? < 1), and Type IV 
_(/3 < -1 but 13 is not an integer)]: T{R,G) = + (3/2) (G/i?). In Ref. [3], 
it has been examine that the finite-time future singularities can be removed 
by the term R^^G^^, where i?i(> 0) and i92(> 0) are positive integers. 

3. Dark energy from F{R) gravity with the Lagrange mul- 
tiphers 

In this section, we study F{R) gravity with the Lagrange multiplier field. 
With Fi{R) and F2{R) arbitrary functions of i?, the action is expressed as 



S = / d^x^ 



Fi{R) - Al ( \d^,Rd^'R + F2{R) 



(10) 



where Al is the Lagrange multiplier field and yields a constraint equation 
(1/2) di^Rdi^R + F2{R) = 0. The variation of the action in Eq. 1^ with 
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respect to g^y leads to the gravitational field equation as 



0, (11) 



where is the covariant derivative and □ = g^'^VfjVu is the covariant 
d'Alembertian. For the de Sitter space-time, which realize the current cos- 
mic accelerated expansion, i.e., the dark energy dominated stage, the scalar 
curvature is a positive constant value Rq and hence the Ricci tensor becomes 
RfMu = (1/4) Rog^u- In this case, from the above constraint equation and 
Eq. (HH), we have Al = [-2Fi(i?o) + Rq {dFi{Ro) / dR)] / [Rq {dF2{Ro) / dR)]. 
Moreover, in the flat FLRW background, the above constraint equation 

reads -{l/2)R? + F2{R) = 0. For F2{R) > 0, this equation can be 
solved in terms of t as t = dR/ y^2F2{R). Provided that the form 
of H{t) is given by the analysis of the observational data, F2{R) is able 
to be reconstructed so that the evolution of H{t) can be reproduced. It 
follows from R = 6 [d/f/dt + 2i?^] that H{t) presents the evolution of 
R = R{t), and by solving this equation inversely, we can find t = t{R). 

Thus, we acquire F2{R) = (1/2) {dR/dtf with t = t{R). We note that 
Fi{R) is an arbitrary function of R. As an example, we consider H{t) = 

ho/t with ho > 1 leading to a{t) = a()t^°, where /iq and ao are con- 
stants. In this case, the accelerated expansion of the universe or power- 
law inflation happens. We have R = 6ho (—1 + 2/io) /t^, from which we 
also acquire t = \/QhQ (—1 + 2ho) /R. Using these relations, we obtain 

F2{R) = R? / [12/io (—1 + 2/io)]. As another example, we examine the case 
that R is described by i? = {R-/2) (1 - tanhwt) + {R+/2) (1 + tanha;t) 
with R±{> 0) and uj(> 0) positive constants. In the limit t — )• ±00, 
R — )■ ±R±, and therefore the universe asymptotically approaches the de 
Sitter space-time. In this case, we can regard that in the limit t — >• 
—00, inflation in the early universe occurs, whereas that in the limit t — t- 
+00, the late-time cosmic acceleration happens. We also have F2{R) = 



(1/8) {R- - R+f u? 1-{R. + R+- 2Rf / {R. - R^ 

quence, for the above R, this F{R) gravity model with the constraint origi- 
nating from the Lagrange multiplier can be a unified scenario between infla- 
tion and dark energy era, although it should carefully be studied whether 
the reheating stage after inflation can be realized. Furthermore, Fi{R) 
does not affect cosmological evolution of the universe and influences only 
the correction of the Newton law. Thus, cosmology is determined only by 
the form of F2{R). 

To explore the Newton law, we take Fi{R) = R/ (2k^) as the Einstein- 
Hilbert term and add matter. In this case, for Al = 0, from Eq. (|11|) 
we find the Einstein equation [Rfj,^ — {1/2) g^^R] = k^T^^ with xjii'P the 
energy-momentum tensor of matter Its trace equation reads R = —k^T'^^\ 



2 



As a conse- 
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where T^'^) is the trace of Moreover, the constraint equation is 

given by (k''/2) d^Td^T + F2 [—k?T^ = 0. Since this is not always met, 
we should modify the constraint equation as (1/2) d^Rd^R + F2{R) — 
(k4/2) 5^r5T(^) - F2 {-k'^T) = 0. Thus, this implies that the action 
with the constraint coming from the Lagrange multiplier field and matter 
should be described by 

S = y^^^^j-^-AL ^d,Rd>^R + F2{R)-'^d^T(^^d>^T('^^ 

-F2(-«2y(^))]+£M}, (12) 

For the case of the vacuum such that T^^^ = 0, the constraint equation is 
(1/2) d^Rdf'R + F2{R) - F2 (0) = 0. If F2 (0) = 0, e.g., the first example of 
F2{R) = R^/ [l2ho (—1 + 2/10)] shown above, this is equivalent to the con- 
straint equation derived from the action in Eq. (jlOp . In this case, there exist 

two types of the solutions in the constraint equation — (1/2) R?'+F2{R) = 0. 

One is -R = and the other is presented by i = dR/ ^J2F2{R). On 
the small scales of, e.g., the solar system and galaxies, the solution would 
be the first solution of i? = so that the Newton law can be recovered. 
On the other hand, in the bulk of the universe, the solution should be 

t = p dR/^j2F2{R) in order that the cosmic evolution can be realized. It 
is not so clear whether the first solution on the small scales of the solar sys- 
tem and galaxies and the second one in the bulk universe can be connected 
in the intermediate scales. 



4. Domain wall solutions in F{R) gravity 

In this section, we investigate a static domain wall solution and reconstruct 
an F{R) gravity model with realizing it [5]. 



4.1. Static domain v^all solution in a scalar field theory 

To begin with, we study a static domain wall solution in a scalar field the- 
ory. We suppose that the following D = d + 1 dimensional warped metric 
ds^ = dy'^ -\-Q^iy) X^^~^o 5^1'^^''^^''' ^^^^ iliat the scalar field only depends 
on y. In this background, the metric in the d-dimensional Einstein man- 
ifold is gpi_y, defined by Rpi_y = \{d — 1) //^] g^y. In addition, for l/f' > 0, 
the space is the de Sitter one, for 1/P > 0, it is the anti-de Sitter, and 
for 1/P = 0, it is the flat. Following the procedure proposed in Ref. [26j, 
it has been demonstrated that a static domain wall solution can exist in a 
scalar field theory [S] (a developed study on a static domain wall solution 
in a scalar field theory has also been executed in Ref. p7j). We inves- 
tigate the action S = J d^ x^f^ [{R/2k^) - {1/2) uJ{Lp)^f,ip^^'ip -V{^)\, 
where uj{ip) is a function of the kinetic term of a scalar field 99 and V((/9) 
is a potential of ip. In the above D dimensional warped metric, with the 
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(y, y) and (/x, u) components of the Einstein equation, we obtain the ex- 
pressions of oj{'^) and V((/5). Using these expressions, the energy density is 

described as p^p = (1/2) uj{ip) {(p') + V((/j). As an example, we consider u = 
uo exp (— y^/yp), where uq and yo are constants. In this case, the distribu- 
tion of reads p^{y) = - [{d - 1) / {2y^)] [(2yVy2) - l] exp (-yV2/o) + 
{d 



1^2 n2 



exp [— no exp (— y^/yp)] . Accordingly, the energy density of 

(fi is localized at y ~ and thus a domain wall is made. We note that a 
condition for p^ to be localized is n — )• in the limit |y| — )• oo. 



4.2. Reconstruction of the form of F{R) 

In the D dimensional warped metric, the (y, y) component and the trace 
of (fj,, v) components of the gravitational field equation read 



d — 1 , , ^ , , d 



Fr -\f = n^T^f) 



d 



P 



fi,iy=0 



(13) 



^'-^e-\Fn 



(14) 



where the prime denotes the derivative with respect to y of d/dy, and 
{F^jlY = dFji/dy and (F/j)" = d^F^ji/dy'^. We examine an explicit form 
of F{R) with leading to a domain wall solution for the case that matter 
is absent. For the model u = exp (— y^/yp) , with the relation R = 

-d^u" + [(1 + d) /4] (u')^} + [did- 1) /P] e"", y can be described as a 

function of R, y = y{R), and eventually we find u = u{y{R)). By plugging 
this equation into Eqs. ([13]) and ([H]) and eliminating y, Eqs. (fT3]) and 
p4p can be expressed as differential equations in terms of F(R). Here, it is 
enough to analyze Eq. (fT3|) because Eq. (fH|) is not independent of Eq. . 
As a result, Eq. (jl3p can be rewritten to 



dR^ 

Ei{R) = {d- 
E2{R) ^ i-d) 



l)u 



,dR 
dy 



dR 

id-1] 



+ 



dR 
dy 



dMyjR)) 

dR^ 



dR 
dy 

= {-d) 

1 fdR 

+ 2 Vdy 



du{yiR)) 
dR ' 
d^Rdu{y{R)) 



dy^ 

2 



dR 



(Mym\ 

\ dR J 



(15) 
(16) 



(17) 
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To solve the above relation of R in terms of y, by defining Y = /uq and 
expanding exponential terms in the limit Y = 'ip' jy^ ^ 1, we take only 
the first leading terms in terms of Y . We find Y = y^/yg ~ (-^ ~ 7i) /72 
with 71 = {2duo/yl) +d{d-l) jl^ and 72 = -d {uo/yl) [6 + (1 + d) uq] + 
\d{d — 1) //^] no, where 71 and 72 are constants. Finally, for Y = y^ ly\ <C 
1, Eq. ^ can be described by {d'^F{R)/dR'^)+C {dF{R) / dR)+VF{R) = 

with C = r!"2^ and V = —\/r!{ \ where 'r!^^ and Hg'^^ constants 
described by the model parameters d, /, uq and yo- We acquire a gen- 
eral solution of this equation as F{R) = Fj^e^+^ + F_e'^~^, where X± = 

(1/2) (^—C ± VC^ — , and F± are arbitrary constants. Here, the sub- 
scriptions lb of \± correspond to the sign "it" on the right-hand side of 
this equation. In the model u = uoexp (— y^/yg), at y ~ the distribution 
of the energy density is localized and therefore a domain wall is realized as 
shown above. Consequently, for an exponential model of F{R) gravity, a 
domain wall can appear at y ~ 0. 

4.3. Effective (gravitational) domain vi^all 

Next, with the reconstruction method |2HI22j. we explore an effective (grav- 
itational) domain wall in F[R) gravity. With the same procedure as in 
Section 2.3, we study the action of F{R) gravity given by F[R, G) = F{R). 
Using two proper functions -P(V') and Q(V') of a scalar field if), we rep- 
resent the term J^{R,G) = P{tp)R + Qitp). The variation over ip yields 
{dP{ip)/dip) R + dQ{ip)/dip = 0. Solving this equation with respect to 
ip leads to ip = ip{R), by substituting which into the action in Eq. ([1]) 
with F{R,G) = P{ip)R + Q{ip), we acquire the action of F{R) gravity as 
F{R) = P{'ip(R))R + Q(^{R)). In the D dimensional warped metric shown 
in Section 4, for the case that ip depends only on y, it follows from the 
gravitational field equation with the choice of ip = y and 1/P = (i.e., the 
flat space), we have 



d-1 



X / dijiPiij))'^^''-^^/^''-'^ {P'{i;)y 



(18) 



d{d-i){u'm 



2 



QW = ^^^P(V) + {d- l)7x'(^)P'(V) , (19) 

where the prime denotes the derivative with respect to V'(= v) of d/dip. 
For a model P{ip) = (C/(V'))"^^'^"^^ and U{ip) = Uq {tp^ + V'o)^ with Uq, ipo 
and X constants, we acquire 

2x^P 32dxV"^-^ 
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t.o f''^ - 1 - 2*) r (2x - 1 - fc) k\ V * 

In the range where = y is large, we take x = ~1/ [4 (4(i — 1)] and impose 
the boundary condition that in the limit \y\ = |^| — >• oo, the universe 
asymptotically approaches flat as u — )• 0. As a result, we obtain u{ip) = 

— {1/ [4 (6d — 1)]} (^o/V') + (^(V'o/V')'^^ • From this expression, we see 

that u{'ip) performs a non-trivial behavior at = y ~ 0. Hence, it can 
be considered that an effective (gravitational) domain wall could appear at 
y = 0. Moreover, by using the representation u'{ip) = {AU' {Tp)/U{ip)) — 
(^Sd/Uiip)"^) J dtp {U' {ip})'^ , we acquire an integration expression of u^ip) as 



ri> I / \ 2{x-i) 

- ^2^^X' / di;— / # + ^2 . (21) 

In Ref. [5j, it has numerically been verified that there exists a local max- 
imum of u{'ip) at V = y ~ 0, and thus an effective (gravitational) domain 
wall could be realized at y = 0. Also, there occurs such a qualitative behav- 
ior of u{'il>) in terms of ip regardless of the values of the model parameters. 

In addition, we mention that for U{^p) = f/o \/V^^~+~0o ' X ~ '^/'^^ there 
exists an analytic solution 




u{ip) = 2(1- 2d) In (V'^ + V'o) + 4d arctan I 1 1 + <^ ' (22 



with C an integration constant. In this case, for the region of a small 
amplitude it is considered that the distribution of the energy density 
is localized, so that an effective (gravitational) domain wall could be made. 

As a demonstration, we reconstruct an explicit F{R) form for u{^lJ) in 
Eq. (j22p . although only in the region of a small amplitude of tp, the distribu- 
tion of the energy density could be regarded as an effective (gravitational) 
domain wall. From P'{il))R + Q'{ip) = and Eq. (fT9|) . we have 

^ = = (d^'(V')n"(^) + 2n"(^)P'(V') + 2n'(^)P"(V')) . 

(23) 

Solving this equation, an analytic relation ip = ip{R) can be found. With 
this relation, we acquire F{R) = P{ip{R))R + Q{'ip{R)). In this case, we 

have P{ip) = {UoipoY^^'^'^^ {l + Yy'^^'^'^\ For Y = ^'^0)1 < 1, by ex- 



panding Eq. (|19p in terms of Y and taking the leading terms, we find R = 
7^o + UiY. Furthermore, from Eq. ([19]) we also obtain Q = QiY + Q2Y'^ ■ 
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Here, TZq, TZi, Qi and Q2 are constants and these are written by using the 
model pararneters d, Uq and j/^q- Moreover, with R = TZq + TZiY we de- 
scribe Y = Yq + YiR, where Yq = —TZq/TZi and Yi = l/TZi. Pi^ip) can also 

be expanded as P{ip) w (C/qV'o)"^^'^"^^ {l - {d - 1)Y + [d{d - I) /2]Y^}. 
By plugging this relation and Q = QiY + Q2Y'^ with Y = Yq + YiR into 
F{R) = P{ip{R))R + Q{iIj{R)) and taking terms of order of i^^ we find 
F{R) = + FiR + where Fq, Fi and F2 are constants and these 

are represented by the model pararneters d, Uq and ifjQ. The above explicit 
form of F{R) has been derived for Y = V'^/V'o ^ ^- Hence, it follows from 
R = TZq + TZiY that this F{R) form can be considered to correspond to the 
one for R ~ 0(1) if TIq ~ 0{1). Therefore, when we choose Fq = {) and 
Fi = 1, we have F{R) = R + F2R^. In this case, for the small curvature 
limit, F(R) approaches R, i.e., general relativity, asymptotically. As a re- 
sult, if u{tP) is given by Eq. (|2'2p in which an effective (gravitational) domain 
wall can be realized, an explicit form of F{R) is expressed as a power-law 
model. We state the difference between the domain walls in Sections 4.2 and 
4.3. A pure gravitational effect yields an effective (gravitational) domain 
wall in Section 4.3, but a scalar field makes a static domain wall solution 
explored in Section 4.1. In Section 4.2, the deviation of F{R) gravity from 
general relativity is equivalent to matter geometrically, i.e., a scalar field in 
Section 4.1. 

5. Covariant higher derivative gravity with scalar projectors 

It is considered that a covariant gravity which is power-counti ng r enormal- 
izable would be higher derivative theory, e.g., models in Ref. [28]. Higher 
derivative gravity is very well known to be renormalizable multiplicatively 
(for a review, see, for example, [29|). However, in general, such a higher 
derivative theory cannot keep the unitarity. To retain it, the so-called 
Hofava gravity |30j has been proposed. In this section, we make the for- 
mulation for covariant higher derivative gravity with Lagrange multiplier 
constraint as well as scalar projectors. In particular, we construct a gravity 
theory with the Lorentz symmetry and/or the full general covariance in 
the action, although these symmetry and/or covariance is spontaneously 
broken. In such a theory, the propagator of the graviton in the ultraviolet 
(UV) region can be improved better, whereas there appears no extra mode 
such as a scalar one. 

5.1. Model 

We explore the following action with the Lagrange multiplier field Al [H |3T] 
as well as the scalar field ^: Sl = - J d^x^f^Xi^ [{1/2) d^^d^"^ + W], 
where ^ and v run 0,1,2,3 and the component denotes the time t as 
Oq = d/dt. From this action, we find a constraint equation (1/2) d^^d^^ + 
W = This means that the vector quantity (5^^*) is time-like one. Hence, 
this breaks the Lorentz symmetry and/or the full general covariance spon- 
taneously. For simplicity, we assume that is a constant, although this 
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assumption is not necessary for the symmetry and/or covariance to sponta- 
neously be broken. Furthermore, the direction of time can be taken so that 
it should be parallel to the vector quantity In this case, we have 

(1/2) (d^/dt)'^ = W, from which we have <I> = V2Wt. Accordingly, the 
spatial region is a hypersurface with a constant <I> because the hypersurface 
becomes orthogonal to the vector quantity (d^^). In the flat space-time, 
we examine the perturbations as = r/^jy -|- with rj^i, the Minkowski 
metric and h^u corresponds to the fluctuations, i.e., the deviation of Qf^i, 
from the Minkowski background r/^,y. A projection operator is defined as 

= ^il' + {df.^d''^) I {2W) with V^^ = 0. As a result, the action de- 
scribing a higher derivative gravity with scalar projector which is covariant 
and power-counting renormalizable is expressed as 



Al ( \d^^d^^ + ) ^ (24) 



52n+3 = y ^'^^^{^ - C [(5'^$5-w^v, - a^ci>5Mw^v,)"p/p, 



Al { \d^,^d^^ + ) ^ (25) 



Here, Eqs. ^ and are for z = 2n + 2 and z = 2n + 3 (n = 0,1,2,. . . ), 
respectively, where z is the quantity denoting the anisotropy between the 
time and spatial coordinates ^30j. Moreover, the gravitational field equa- 

tion is given by [1/ (2^2)] [R^, - (1/2) g^.i?] + 4'^''''''"^ " (Al/2) + 
(1/2) g^^ [(1/2) dp<^>dP<^> + W]=0 with G^ftu^^'''^ the higher derivative term, 
i.e., the second term, in the actions in Eq. (|24]l and ([25]) . Suppose the flat 
vacuum solution, the constraint equation (1/2) d^^^d'^^ + W = becomes 

(1/2) {d^/dty = W. For the flat space solution, the gravitational field 
equation is reduced to Xi,dfj,^du^ = 0, because all of the term V^^V^^ 
as well as the curvature terms vanish. The solution is given by Al = 0, 
because 9^$ / owing to (1/2) d^^^d'^^ + W = 0. Thus, in these actions 
in Eqs. ()24p and (j25p solutions with Al = in the flat space vacuum can 
be realized. We further analyze the perturbations hfj_u with Al = 0. Using 
the diffeomorphism invariance in terms of the time coordinate, as a gauge 

condition for the unitarity, we set <I> = V2W t. By taking only the quadratic 
terms of the perturbations, we rewrite the actions in Eqs. (j24p and (j25p . 
At this stage, there remains the diffeomorphism invariance in terms of the 
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spatial coordinates. In addition, the term hoi in the higher derivative term 
with a coefficient does not exist in the rewritten actions. Moreover, the 
above constraint equation leads to /iqo = 0. We derive the equations by 
varying the actions in Eqs. (pH) and ([25]) with respect to /iqo and ^. We 
decompose /loi as h^idiS + Vi with d^Vi = 0, where s is the spatial scalar 
quantity and Vi is a vector field. The invariance in terms of the spatial 
coordinates under the transformations of the linearized diffeomorphism is 
described as 6x^ = d^u + with diW^ = 0, where u is the spatial scalar 
quantity and Wi is a vector field. We find the following transformations un- 
der the diffeomorphism: 5s = dtu and 6vi = dfWi. Accordingly, the gauge 
condition s = f* = 0, i.e., hu = can be chosen. Furthermore, we express 
hij as hij = 5ijA + djBi + diBj + Cij + [did j — (1/3) 6ijdkd^] E, where A 
and E are scalar quantities, Bi is a vector field, and Cij is a tensor field. 
Here, d'Bi = 0, d'Q 



0, and Cj* = 0. As a result, we acquire 
Al = and vector Bi = 0, and thus the scalar Al and vector Bi modes do 
not propagate. We fix the gauge in the actions rewritten above. We vary 
these actions over A and E and obtain equations. Prom these equations and 
A = (1/3) dkd^E, which is derived by the equation derived by the variation 
of the action with respect to /iQi with the above decomposed representation 
of hij, we find OqA = 0. Hence, since A and E have only the dependence 
on the spatial coordinate, these scalar quantities do not propagate. As a 
consequence, all the scalar modes Al, /iqO) s, A and E and all the vec- 
tor modes Vi and Bi do not propagate, whereas the propagating mode is 
only the tensor mode Cij, namely, a massless graviton. This is a different 
feature from the Hofava gravity [3D] without the Lorentz invariance. The 
final expressions of the actions in Eqs. (j24p and ()25p are given by 



1 

8^ 



( -di + dkd" ) c'^ 



n+1 



n+l 



S- 



2n+3 



L 



C^J ( -dl + dkd'' ] 



)2n— 1 ATi/2n+l 



CIV' 



n+l 



C,;. 



n+2 



(26) 



(27) 



Therefore, in the momentum space the propagator reads 



{hij{p)hki{-p)) 



PiPj 

p2 



PkPi^ _ 



- 

1/ 
1/ 



PiPi 

p2 



5, 



p. 

PjPk \ 

p2 ; 



PiPk 

p2 



hi 



PjPi] 



2 j 



22'^CK;2^2np4(n+l)j 
22n-l (A^2 y^2n+lp2{2n+3) 1 



for 
for 



2n + 2 
2n + 3 



(28) 
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with = Y^i=i {P^Y ^^"^ = — {p^Y + P^- If C > and = 0, when 

22"(K2Ty2np4n+2 = 1 for z = 2n + 2 and 22«-l(K2^2n+lp4{n+l) ^ ^ 

z = 2n + 3 are satisfied, the tachyonic pole exists. Thus, at least the flat 
vacuum is unstable. In this model, at least on the tree level, no propagating 
vector or scalar mode exists. The fact that the tensor structure of the 
propagator in Eq. (j28p changes implies that the vector or scalar mode could 
emerge. In other words, the vector or scalar mode has to be a composite 
state. At any perturbative level, this does not appear usually. Accordingly, 
the quantum corrections should not change the tensor structure. In the UV 
region with large fc, for z = 2 (n = 0) in Eq. (j26p . the propagator evolves 

as l/|fc| . Hence, the UV behavior performs. While, for 2: = 3 (n = 0) 

in Eq. ([27|) . the propagator evolves as 1/ |fc| . Thus, the model is power- 
counting renormalizable. For 2; = 2n + 2 (n ^ 1) in Eq. ()26p or z = 2n + 3 
(n ^ 1) in Eq. (|27p . the model is power-counting super-renormalizable. 
In the high energy regime, the dispersion relation of the graviton becomes 
uj = ck^ with c(> 0) a positive constant for the consistency of the dispersion 
relation. Here, u is the angular frequency which corresponds to the energy 
and k is the wave number which does to the momentum. Accelerating 
cosmology in such a theory was studied in Ref. [6] . 

6. Conclusions 

We have studied the accelerating (dark energy) solutions of modified grav- 
ity. These solutions may yield future singularities. We have explored the 
finite-time future singularities in F{R), F{G) and J-{R, G) gravity theories. 
The removal of the finite-time future singularities in F{R) gravity by adding 
an E?-ieim which simultaneously leads to the unification of early-time infla- 
tion with late-time acceleration [9] has been mentioned. The corresponding 
term may be different for F{G) or J^(R,G) gravity theory |3j. Moreover, 
we have studied dark energy in F{R) gravity with the Lagrange multi- 
plier field. Furthermore, domain wall solutions in F{R) gravity have been 
presented. In addition, we have investigated covariant higher derivative 
gravity with scalar projectors. 
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